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EFFECTS  OF  ION  COLLISIONS  ON  QUASILINE AR  HEATING 
BY  THE  CURRENT-DRIVEN  ION  CYCLOTRON  INSTABILITY 
IN  THE  HIGH  LATITUDE  IONOSPHERE 

I.  INTRODUCTION 

Recently,  much  interest  has  been  devoted  to  the  role  of  the  high 
latitude  ionosphere  in  influencing  oagnetospheric  dynamics  and  morphology. 
In  particular,  the  problem  of  the  high  latitude  ionosphere  as  a  source  of 
magnetospheric  plasma  has  come  under  intense  study  in  the  last  few  years. 
Several  investigators  have  observed  upvard  flowing  accelerated  ionospheric 
ions  in  the  auroral  zone  space  plasma  [Shelley  et  al.,  1976;  Sharp  et  al., 
1977;  Mizera  and  Fennel,  1977;  Ghielmetti  et  al.,  1978;  Croley  et  al., 
1978;  Whalen  et  al.,  1978;  Klumpar,  1979;  Collin  et  al.,  1981;  Gorney  et 
al.,  1981;  Yau  et  al.,  1983;  Hoore  et  al.,  1986].  The  particle 
distributions  of  these  upf loving  ions  can  be  characterized  as  either  beam¬ 
like  with  pitch  angles  near  180  degrees  or  conical  with  pitch  angles 
betveen  90  and  approximately  130  degrees.  Gorney  et  al.  [1981]  have 
concluded  that  the  beam-like  distributions  are  most  often  observed  above 
4000-5000  km.  Conic  distributions,  on  the  other  hand,  have  been  observed 
near  1500  km  [Klumpar,  1979;  Ungstrup  et  al.,  1979]  and  at  rocket  altitudes 
near  400-600  km  [Whalen  et  al.,  1978;  Yau  et  al.,  1983,  Basu  et  al.,  1988] 
suggesting  a  low  altitude  generation  region  for  ion  conics.  Furthermore, 
Gorney  et  al.  [1981]  have  shovn  that  the  large  majority  of  conic 
distributions  have  energies  below  a  few  hundred  eV  implying  a  significant 
ionospheric  component.  Moore  et  al.  [1986]  have  recently  observed 
perpendicular  ion  heating  at  high  (700  km)  and  low  (300-350  km)  altitudes 
with  the  TOPAZ  sounding  rockets.  Their  observations  show  that  ion  heating 
at  700  km  is  correlated  with  the  lower-hybrid  wave  activity  while  ion 
heating  observations  at  300  km  is  correlated  with  keV  electron 
precipitation.  They  do  not  observe  any  significant  ion  heating  below  270 
km. 
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Kintner  et  al.  [1978],  using  S3-3  satellite  data,  have  measured 
electrostatic  ion  cyclotron  waves  near  regions  of  upward  flowing  ions.  At 
lower  altitudes  (400-600  km)  Yau  et  al.  [1983]  have  observed  ion  cyclotron 
wave  activity  near  ion  transverse  acceleration  regions.  Bering  [1984]  has 
given  evidence  of  ion  cyclotron  wave  activity  also  at  low  altitudes  (350 
km)  in  the  diffuse  aurora.  Recently,  Basu  et  al.  [1988]  have  observed 
apparent  ion  conics  at  approximately  350  km  in  the  presence  of  wave 
activity  in  the  ion  cyclotron  range  of  frequencies  and  large  field  aligned 
currents.  In  addition  to  low  altitude  ion  heating  in  the  ionosphere, 
intense  local  heating  has  been  observed  at  much  higher  altitudes  in  the 
magnetosphere  where  the  plasma  is  collisionless.  Recent  papers  have  also 
proposed  that  ion  heating  could  occur  cumulatively  over  several  thousand 
kms  [Temerin,  1986;  Retterer  et  al.,  1987;  Andre  et  al.,  1988]. 

The  ion  conic  distributions  have  been  interpreted  as  resulting  from 
perpendicular  heating  of  ions  at  low  altitudes  followed  by  parallel  upward 
adiabatic  motion  due  to  the  magnetic  mirror  force  [Sharp  et  al.,  1977]. 
Proposed  perpendicular  heating  mechanisms  include  collisionless 
electrostatic  ion  cyclotron  waves  [Ungstrup,  1979;  Okuda  and  Ashour- 
Abdalla,  1981],  lower  hybrid  waves  [Chang  and  Coppi,  1981],  and 
electrostatic  shocks  and/or  oblique  double  layers  [Mozer  et  al.,  1980;  Yang 
and  Kan,  1983].  Recently,  Satyanarayana  et  al.  [1985]  have  discussed  the 
linear  theory  of  the  electrostatic  ion  cyclotron  instability  in  the  high 
latitude  ionosphere  where  collisional  effects  are  important.  They  have 
shown  that  electron  collisions  are  destabilizing,  under  realistic 
ionospheric  conditions,  and  lead  to  the  excitation  of  the  collisional  ion 
cyclotron  instability  [D'Angelo,  1973;  Chaturvedi  and  Kaw,  1975; 
Chaturvedi,  1976;  Suszcynsky  et  al.,  1986]  in  the  bottomside  ionosphere. 
Ion  collisions  were  found  to  have  a  stabilizing  influence  on  the 
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collisions!,  ion  cyclotron  instability.  Dupree's  resonance  broadening 
theory  was  used  by  Satyanarayana  et  al.  [1985]  to  estimate  the  nonlinear 
saturated  aaplitudes  of  the  collisional  ion  cyclotron  instability  in  the 
high  latitude  ionosphere.  Furthermore,  Satyanarayana  et  al.  [1985] 
conjectured  that  the  collisional  ion  cyclotron  instability  should  be  a 
source  of  transversely  accelerated  heavy  ions  at  higher  altitudes  as 
observed  in  several  studies  [Hoffman  et  al.,  1974;  Hultqvist,  1983;  Craven 
et  al.,  1985]. 

Hovever,  Satyanarayana  et  al.  [1985]  did  not  discuss  the  ion  heating  by 
ion  cyclotron  vaves  in  regions  of  the  collisional  high  latitude  ionosphere 
nor  did  they  discuss  the  regions  of  the  high  latitude  ionosphere  where  ion 
conic  formation  might  be  possible.  In  this  brief  report  we  study  the 
effects  of  ion-ion  and  ion-neutral  collisions  such  as  isotropization  and 
cooling  on  quasilinear  ion  heating  by  the  current-driven  ion  cyclotron 
instability  in  the  high  latitude  ionosphere.  Ve  model  the  heating  with 
collisionless  quasilinear  heating  rates  to  determine  the  threshold  altitude 
for  conic  observations  in  the  bottomside  ionosphere.  In  Section  2  we  give 
the  basic  collisional  model  used  for  the  ions  in  this  study.  In  Section  3 
we  discuss  approximate  quasilinear  heading  rates  with  collisional  effects. 
Finally,  in  Section  4  we  summarize  and  discuss  our  results. 

II.  MODEL 

tfe  consider  a  homogeneous  collisional  plasma  in  a  uniform  magnetic 

field  B  along  the  z-direction,  i.e.,  B  =  Bqz.  In  this  study,  ve  use  local 

theory  and  ignore  inhomogeneous  effects,  e.g.,  density  and  temperature 

gradients.  As  a  result,  we  assume  k^L^  >>  1  and  k^L^  >>  1  where  and  Ln 

a  L  are  the  gradient  scale  lengths  of  temp-  ature  and/or  density  in  the 
z 

direction  perpendicular  and  parallel  to  the  magnetic  field,  respectively, 
and  k^  and  k()  *  kz  are  the  perpendicular  and  parallel  components  of  the 
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instability  vavevector  k.  We  consider  the  altitude  range  of  approximately 
200-600  ka  in  the  high  latitude  auroral  and  polar  ionosphere.  In  this 
altitude  regime  electron  and  ion  collisional  effects  are  important 
[Satyanarayana  et  al.,  1985 J  for  the  evolution  of  the  current-driven  ion 
cyclotron  instability.  The  electrons  are  treated  in  the  collisionless 

liait  in  this  study.  We  consider  only  ion-ion  and  ion-neutral  collisions 

1/2 

since  for  200  <  z  <  1000  ka  \^n,  >  \>^e  since  \>„  =  (m^/m  )  \c  , 
[Spitzer,  1962;  Trubnikov,  1965],  where  is  the  ion-ion  collision 
frequency,  \^n  the  ion-neutral  collisional  frequency,  the  ion-electron 
collision  frequency,  the  ion  mass  and  m^  the  electron  mass. 

The  basic  equation  describing  the  evolution  of  the  plasma  distribution 
function  f  of  the  current-driven  electrostatic  ion  cyclotron  instability 
in  a  collisional  plasma  can  be  written  as 


3f  e  ,  ,  .  ,3f  . 

t-2  +  v  •  7f  ♦  — |E  ♦  ±v  x  B|  •  7  f  = 

3t  -a  a  m  i.-  C-n  - J  v  a  ^3t  )( 


(1) 


where  a  denotes  species  i  or  e,  e  the  charge  of  species  a,  m  the  mass  of 

a  at 

species  a,  E  is  the  electric  field,  B  is  the  magnetic  field,  c  the  speed 

of  light,  7  ■  3/ 3v  ,  7  ■  3/3x,  and  (3f  /3t)  is  the  contribution  due  to 

collisional~?ffects.  Satyanarayana  et  al.  [1985]  has  given  a  complete 

discussion  of  the  linear  theory  of  the  current-driven  electrostatic  ion 

cyclotron  instability,  using  Eq.  (1),  with  (3fa/3t)c  given  by  a  BGK 

collision  operator  both  for  electrons  and  ions.  In  this  study  we  take  a 

Fokker-Planck  representation  [Montgomery  and  Tidman,  1964]  of  (3fa/3t)c  for 

ion-ion  Coulomb  collisional  effects.  This  representation  can  be 

approximated  by  using  an  effective  relaxation  time  model  [Trubnikov,  1965]. 

Por  ion-neutral  collisions  we  adopt  the  approach  outlined  in  Satyanarayana 

et  al.  [1985]. 
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III.  QUASILINEAR  ION  HEATING  WITH  COLLISIONAL  COOLING 

Ve  take  the  ion  distribution  function  to  be  a  bi-Maxvellian  with 
effective  teaperatures  T^  and  T|(  ■  Tz  in  the  direction  perpendicular  and 
parallel  to  the  magnetic  field: 


.  (  )  f  "i  u  i  n/2  r  vi  vu 
filVWllJ  ■  UnTnJ  exp  [~  20^/10.)  '  2(TM/m.) 

In  lieu  of  an  exact  treatment  of  quasi-linear  ion  heating  including  ion-ion 
collisions,  as  implied  by  (1),  ve  make  an  approximation.  In  the  limit  of 
weak  collisions  ve  set  the  right  hand  side  of  (1)  to  zero  and  solve  for  the 
quasi-linear  ion  heating  rates  vithout  ion  collisional  effects  (first  term 
on  right  hand  side  of  (2)  and  (3)  below)  as  has  been  derived  by  others 
[Davidson,  1972;  Dakin  et  al.,  1976.,  Ashour-Abdalla  and  Thorne,  1978, 
Okuda  and  Ashour-Abdalla,  1981].  In  the  next  order  ve  neglect  the  wave 
effects  and  solve  (1)  for  the  effects  of  ion  collisions.  The  principal 
ion-ion  collisional  effect  gives  a  temperature  isotropization  term 
proportional  to  (T^  -  T|()  using  an  effective  relaxation  time  model 
[Trubnikov,  1965]  (second  term  on  right  hand  side  of  (2)  and  (3)).  The 
quasilinear  heating  rates  vith  ion-ion  collisions  in  a  single  species 
plasma  can  then  be  written: 


x  exp 


[- 


n2.) 


2k^T 


''efflv  Tn)(Ti  -  T„) 


(2) 
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(3) 


dTll 

dT 


„  n<“  *  n2i>rn(s)l\1' 
^  k,, 


x  «P  I"  JL  ]  +  2veff ^Tl’  Tl,><Tl  -  V 

ZHI  Cll  ’ 

2  2 

vhara  Qj  -  eflg/BjC,  E  «  -  7+,  Tn(s)  »  Tn(s)exp(-s) f  s  ,  k^,  the 

modified  Bessel  function  of  order  n,  is  the  linear  frequency, 

Fourier  transform  of  the  electrostatic  potential  t  such  that  e^/Te  is 

dimensionless,  v#ff  -  2n1/2e4n.Xm.'1/2(T||)"3/2A'2[ (A  ♦  3)  ( tan_1A1/2)/A1/2 

-  3]sec  A,  A  ■  23  -  ln((n/T;A  ]  is  the  Coulomb  logarithm  and  A  *  (T^  - 

T||)/T||  >  0.  tfe  note  that  for  »  0  the  quasilinear  heating  rates  reduce 

ta  those  given  in  Eq.  (3)  and  (6)  in  Okuda  and  Ashour-Abdalla  [1981]. 

Equations  (2)  and  (3)  describe  the  evolution  of  T^(t)  and  T|(( t )  with 

quasilinear  heating  as  a  source  and  ion  collisional  cooling  as  a  sink. 

(Eq.  (2)  and  (3)  include  only  the  effects  of  ion-ion  collisions  in  the  form 

of  In  Eq.  (6)  and  (7)  ve  generalize  the  model  and  add  a  second  ion 

species  and  ion-neutral  collisions  as  veil.)  The  approximate  quasi-steady 

state  ion  temperatures  and  T|(  and  anisotropy  (T^  -  Tu)  can  be  found  by 

balancing  the  source  and  sink.  An  analysis  similar  to  ours  has  been 

performed  by  Ionson  et  al.  (1976)  who  considered  collisional  effects  on  ion 

cyclotron  wave  heating  at  altitudes  greater  than  2000  km  in  the  auroral 

zone.  However,  they  did  not  compute  a  minimum  altitude  below  which  ion 

cyclotron  wave  heating  would  be  isotropized  by  collisions  nor  did  they  use 

quasilinear  effects  as  an  ion  heating  source. 

Equations  (2)  and  (3)  are  difficult  to  solve  exactly  since  the  spectrum 
2 

1^1  of  the  current-driven  ion  cyclotron  instability  is  not  well  known. 


is  the 
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Numerical  simulations  (Okuda  and  Ashour-Abdalla,  1981]  indicate  that  |<^| 
peaked  about  k  corresponding  to  the  fastest-groving  fundamental  «  =  2^ 
harmonic  mode.  Equations  (2)  and  (3)  can  be  written,  approximately, 


3Ti  - 

fr%2 

!s!fl2  ri(s)f  "i  H/2 

f  -ii 

3 r " ' 

l  Te  J 

k(|  l2nTnJ  exp 

L  2kf|T,,  Jj 

-  WV  T>I)  (T1  -  V  (4) 


|JS 

1 

(  e\ 

e 

1  Fl<s> 

(  rai  11/2 

exp 

f  "l<\  ‘ ,1 

at 

1  k„ 

+  2\>  ,, 
eff 

(Tr  T„ 

)  (T±  -  T,,) 

(5) 

where  k  corresponds  to  the  instability  wavevector  k  giving  maximum  linear 
growth  for  the  kinetic  collisionless  ion-cyclotron  instability.  In  order 
to  model  a  more  realistic  ionosphere,  we  generalize  (4)  and  (5)  by  adding  a 
second  minor  species  and  ion-neutral  collisions  [Satyanarayana  et  al., 
1985].  The  ion-neutral  collisions  help  model  the  eventual  wave  energy 
transfer  to  neutrals.  We  model  the  multispecies  collisional  effects  with  a 
collision  frequency  [Banks  and  Kockarts,  1973]  and  ion-neutral 

collisions  with  the  collision  frequency  \k  .  In  a  multicomponent  plasma, 
e.g.,  N0+  and  0+,  Eq.  (4)  and  (5)  become 

TT  -  °<Ti'  Vai  -  Vf<Ti-  V<Ti  -  tn>  -  VT1  -  V  -  -  V 

(6) 
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3T, 


if  “  Q<Ti’Tn><«k  -  ai>  +  2^eff<T±*Tn><Ti-Tn>  -  VW  -  "m<  W 


(7) 


where 


0(Ti’  V  -  K  /Tj2  <Te/Ti>ri<sl>[ai("i/2',rii,1''2/kli]e’"'[--i57^  ] 

-m  ikuTn 


/  .1/2  , 

(m  m  )  n  X 

.  „  . «-19  N0+  0+  0+  N0+0+  -1 

v  >  1.8  x  10  - ra-rs—  sec 

m  L  m  _  „  U/2 


[m  T  +  m  T  S' 
^  N0+  0+  0+  N0+^ 


-21  1/2  -1 
v.  »  6  x  10  a  n  (T./m.)  sec  1 
in  n  l  l 


0^  is  the  ion  cyclotron  frequency  of  the  majority  ions,  nn  is  the  neutral 
density,  is  the  temperature  of  minority  ions  (^Nq+,Tq+  or  TH+)  assumed 


H 


to  be  in  equilibrium  with  the  electrons  (T^  -  Tg),  A  is  the  Coulomb 

logarithm,  Tn  the  neutral  temperature,  and  T|(  refer  to  the  temperature 

perpendicular  and  parallel  to  the  geomagnetic  field  of  the  majority  ion 

species,  with  all  other  symbols  retaining  their  standard  meaning. 

In  the  following,  we  fix  (e^  /Te>,  kt|/k^  and  solve  (6)  and  (7)  for 

-in  ~ 

T^(t),  T||(t).  For  simplicity  of  notation,  in  the  following  we  drop  the 

subscript  k^.  We  estimate  typical  values  for  e$/Te>  From  Yau  et  al  [1983] 

and  Bering  [1984]  we  use  typical  electric  fields  E  =  10-40  mV/m  in  regions 

of  ion  energization  and  find  e$/T  =  eE/k  T  =  0.2  -  0.4  for  T  =  3000°K 

erne  e 

and  k^  a  p~*  with  p.  =  2.5m  the  approximate  0+  ion  gyroradius  in  the  high 

latitude  ionosphere.  This  range  of  values  of  e$/Te  is  also  consistent  with 
computer  simulations  [Okuda  and  Ashour-Abdalla,  1981;  Pritchett  et  7., 
1981]  of  the  current  driven  electrostatic  ion  cyclotron  instability  and  S3- 
3  satellite  observations  at  high  altitudes  [Mozer  et  al.,  1980;  Kintne.  et 
al.  1979].  The  quantity  (w  -  2^)  is  taken  to  be  -  2^  from  the  linear 
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dispersion  relation  of  current  driven  ion  cyclotron  waves  (Satyanarayana, 
et  al. ,  1985).  We  set  k|(  *  0.1  k^  since  this  value  of  k.()  gives  maximum 
linear  grovth  [Satyanarayana  et  al. ,  1985].  In  addition  ve  normalize  T^, 
T|(  vith  respect  to  Tg  and  the  time  t  is  normalized  by  2^  where  2^  is  the 
ion  gyrofrequency  of  the  majority  ion. 

Pig.  1  gives  the  evolution  of  and  T(|  according  to  Eq.  (6)  and  (7)  as 
a  function  of  time  at  several  different  altitudes  in  the  high  latitude 
ionosphere  given  in  Table  1  [Banks  and  Kockarts,  1973,  Schunk  et  al.,  1975, 
1976].  In  Pig.  1  we  fix  e$/Tg  -  0.2.  Curves  A,  B,  C,  and  D  correspond  to 
altitudes  200,  300,  4C0,  and  600  km,  respectively.  In  general,  ve  find 
increasing  anisotropy  with  altitude  vith  T^/T^  =  1.4  at  low  altitudes 
(curve  A,  200  km)  and  T^/T^  =  3.1  at  higher  altitudes  (curve  D,  600  km)  for 
several  hundred  ion  cyclotron  periods. 

Pig.  2  displays  the  evolution  of  T^(t)  and  T|((t)  at  the  same  altitudes 
as  Fig.  1  but  with  a  larger  ion  cyclotron  wave  amplitude  e4/Tg  =•  0.4. 
Curves  A,  B,  C,  and  D  correspond  to  the  same  altitudes  as  in  Fig.  1.  From 
Pig.  2  we  find  increased  heating  compared  to  the  smaller  wave  amplitude 
case  in  Pig.  1.  We  also  note  increasing  anisotropy  with  altitude  with 
T^/T||  a  1.9  at  low  altitudes  (curve  A,  200  km)  and  T^/^  =  6.4  at  higher 
altitudes  (curve  D,  600  km)  again  on  the  time  scale  of  several  hundred  ion 
cyclotron  periods.  The  evolution  of  T^(t)  and  T|((t)  using  variations  of 
the  model  ionospheric  parameters  contained  in  Table  1  was  also  studied  vith 
the  results  qualitatively  similar  to  those  displayed  in  Figs.  1-2. 

The  threshold  altitudes  above  which  ve  find  strong  anisotropic  heating 
are  near  or  just  below  the  F-region  peak  in  the  high  latitude  ionosphere 
for  wave  amplitudes  a$/Tg  =  0.2,  0.4,  respectively.  We  note  that  both  the 
ion-neutral  and  ion-ion  collisions  contribute  to  the  isotropization  below 
the  F-region  peak.  Since  heavy  ions  such  as  NO^O^)  or  0+  are  major 
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species  at  these  altitudes,  ve  conjecture  that  heating  by  large  amplitude 
collisional  EIC  may  cause  heavy  ion  conic  distribution  above  the  F  region. 


IV.  SUMMARY  AND  CONCLUSIONS 

Using  an  effective  relaxation  time  model,  ve  have  made  a  preliminary 

study  of  the  effects  of  ion-ion  and  ion-neutral  collisions  on  quasilinear 

ion  heating  due  to  the  current-driven  ion  cyclotron  instability  in  the  high 

latitude  ionosphere.  For  conditions  typical  of  the  upper  altitude  high 

latitude  F-region  ionosphere  (600  km)  ve  find  that  strong  temperature 

anisotropies  T^/T^  >3  for  ')+  can  be  sustained  even  in  the  presence  of  ion 

collisional  effects  for  several  hundred  0+  ion  cyclotron  periods  for  ion 

cyclotron  vave  amplitudes  e$/Tg  -  0.2  and  0.4  as  predicted  by  quasilinear 

theory.  For  the  lover  altitude  F-region  ionosphere  (200-300  km)  ve  find 

approximate  collisional  isotropization  T^/T( (  a  1  for  vave  amplitudes  et/T^ 

■  0.2  on  time  scales  of  a  fev  tens  to  hundreds  of  N0+  ion  cyclotron 

periods.  For  larger  vave  amplitudes  e+/Te  *  0.4  ve  find  only  veakly 

anisotropic  heating  at  lover  F-region  altitudes  (200-300  km).  We  conclude 

that  ion  conics  driven  by  the  electrostatic  ion  cyclotron  instability  vill 

be  seen  at  altitudes  greater  than  approximately  300-400  km  in  the  high 

latitude  ionosphere.  This  finding  is  consistent  vith  recent  data  derived 

from  rocket  experiments  [Moore  et  al.,  1986]  vhich  indicate  significant  ion 

perpendicular  heating  at  altitudes  above  250-300  km  in  the  high  latitude 

ionosphere.  Our  results  are  also  consistent  vith  recent  DE  satellite 

observations  of  Basu  et  al.  [1988]  vho  measure  conic  formations  only  for 

5  -3 

densities  n@  <  2  x  10  cm  in  the  presence  of  large  field-aligned  currents 
and  ion  cyclotron  vave  activity.  In  addition,  it  should  be  pointed  out 
that  Figs.  1  and  2  shov  perpendicular  ion  heating,  even  in  those  cases 
vhere  near  isotropy  occurs.  It  is  clear  from  Figs.  1  and  2  that  curves  B-D 
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(altitudes  from  300-600  km)  show  perpendicular  ion  heating  factors  from  2- 
10  depending  on  altitude  and  EIC  wave  amplitudes  (larger  perpendicular 
heating  occurring  at  higher  altitudes  and  larger  wave  amplitudes).  For 
Fig.  2  (e+/Tg  -  0.4)  even  curve  A  (altitude  of  200  km)  shows  a 
perpendicular  ion  heating  factor  -2.  Consequently,  if  we  assume  the 
quasilinear  heating  rate  to  apply  for  the  collisional  EIC,  collisional  ion 
cyclotron  waves  could  provide  a  means  for  a  low  altitude  pre-acceleration 
(or  heating)  region. 

In  this  study  we  have  used  am  essentially  empirical  model,  the 
effective  relaxation  time  model,  to  study  the  effects  of  ion  collisions  on 
quasi-linear  heating  by  the  current-driven  ion  cyclotron  instability  in  the 
high  latitude  near  earth  space  plasma  environment.  Our  aim  was  to  develop 
a  simple  model  of  time  scales  for  ion  collisional  relaxation  in  the  high 
latitude  ionospheric  environment  and  to  balance  these  relaxation  times  with 
those  associated  with  wave  heating  in  order  to  better  characterize  possible 
ion  conic  generation  and  formation.  In  order  to  gain  more  complete 
information  about  the  ion  distribution  function,  collision  models  that  are 
both  more  realistic  and  sophisticated  will  have  to  be  employed.  We  defer 
such  studies  to  future  work. 
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Table  1 


Curve 

Altitude(km) 

n^Ua"3) 

_3 

n0+(cm  ) 

_3 

nn(cm 

A 

200 

1  x  105 

1.5  x  104 

7.5  x  10 

B 

300 

2.1  x  104 

1  x  105 

7.5  x  10: 

C 

400 

2.5  x  103 

7  x  104 

1.2  x  101 

0 

600 

1  x  102 

2.8  x  104 

5.6  x  10 

Note:  ne  -  ♦  V 
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